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We present our methods to fit the two point correlators for light, strange, and charmed pseu- 
doscalar meson physics with the highly improved staggered quark (HISQ) action. We make use 
of the least-squares fit including the full covariance matrix of the correlators and including Gaus- 
sian constraints on some parameters. We fit the correlators on a variety of the HISQ ensembles. 
The lattice spacing ranges from 0.15 fm down to 0.06 fm. The light sea quark mass ranges from 
0.2 times the strange quark mass down to the physical light quark mass. The HISQ ensembles 
also include lattices with different volumes and with unphysical values of the strange quark mass. 
We use the results from this work to obtain our preliminary results of fo, /b s , fo JfD, and ratios 
^ of quark masses presented in another talk [1]. 
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1. Introduction 



Fitting correlators is one of the essential procedures in studying non-perturbative physics in 
Lattice QCD. The hadron spectrum and decay constants can be calculated by fitting two-point 
correlation functions, and hadronic form factor and matrix elements can be calculated by fitting 
three -point functions. Here we present our method for fitting two-point meson correlators with 
HISQ. 

The meson correlators are computed using the HISQ action[4], which is an &{a 2 ) -improved 
staggered quark action, with smaller taste-symmetry breakings than in the asqtad action we used 
previously. For the pseudoscalar meson correlators, we use both random-wall sources and Coulomb 
gauge fixed wall sources, with the point operator as the sink. We do combined fits to the correlators 
of both sources. This helps us to isolate the ground states that we are interested in. More details 
about the correlators are given in Sec. 2. 

The HISQ lattices include the vacuum polarization of four dynamical quarks: up, down, 
strange, and charm quarks [2, 3]. There are, in total, nineteen HISQ ensembles with different 
lattice spacings, quark masses, and spatial volumes. The parameters of the HISQ lattices are tabu- 
lated in Table 1 . Note that we have physical sea-quark mass ensembles, where all sea-quark masses 
are tuned to take the physical values. These ensembles eliminate the need for a chiral extrapola- 
tion, and they can be used, together with higher unphysical-mass ensembles, to correct for slight 
mistunings of the sea quark masses. 

Since elements of the correlators at different time slices are highly correlated, we use the full 
covariance matrix, which includes both diagonal and off-diagonal elements, to define the objective 
function % 2 . We use Gaussian constraints, or priors in the language of Bayesian fits, on some 
parameters of heavier states. This is useful especially when multiple states are needed to get reliable 
fits. 

We divide the correlators into two groups: light-light and heavy-light. In this terminology, 
light refers to the up, down, and strange quarks, while heavy refers to the charm quark. We find that 
we can get good fits with a single state for light-light correlators as long as large enough distances 
are chosen for the fitting range. However, we find multiple state fits are efficient for heavy-light 
correlators to get good fits with small statistical errors. Our fitting methods are discussed in more 
detail in Sec. 3. 

We block the correlators in the Monte Carlo trajectory to account for the autocorrelations. To 
find an optimal size of the blocks, we look at how the covariance matrix scales as the size of the 
blocks is increased. In the course of this analysis, we find that the autocorrelations depend on 
(Euclidean) time, and that they can be understood as a consequence of the way the source time 
slices are chosen. A short discussion about this is given in Sec. 4 

The systematic error due to the excited-state contamination is tested by varying the fitting 
ranges and priors. Changes from these variations are small compared to our statistical errors. This 
is briefly discussed in Sec. 5 
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2. Meson Correlators 

We use the point and wall operators to create the pseudoscalar mesons, 

Op(x,t) = ZA(x,t)e(x,t)zB(x,t) , (2.1) 
Ow(f) = J^ZA(x,t)e(x + y,t)xB(y,t) , (2.2) 

x.y 

where % is the HISQ field, and A, B are flavor indices. e(x,t) = (— iy^i x < +t corresponds to the 
projection to the Goldstone pion with staggered quarks [5]. The Coulomb gauge fixing must be 
used for the wall operator, so we call this Coulomb wall operator. 

Using the point sink, we can define two correlators depending on the source, 

Cw{t) = (yJ^O P (x,t)ol{0)) . (2.4) 

The summation over the sources for the point operator is implemented using random noise, so we 
call Cr random wall source correlators. 

These correlators can be expressed in terms of masses and amplitudes of relevant excitations. 
We adopt the following parameterization. 

C R (t) = YAM] ( e-^' + e-^-'A + K Y(-\) t A' k M'l ( +e -M' k (Nr-t)\ ^ (2 5) 
C w (f)= ^BjMUe-^' + e-^-'A + £ (-1)'*^ (e'^ + , (2.6) 

j=0 ^ ' k=Q ^ ' 

where / and K are the numbers of ordinary and alternate states respectively. The backward prop- 
agations from the image source at Nj, the temporal size of the lattice, are included because of 
periodic boundary condition on the mesons. 

Following [7, 8], we fit the two correlators simultaneously with common masses. It helps us to 
isolate the ground states of the point source since the Coulomb wall operator is less contaminated 
from excited states. Once the masses and amplitudes are fitted, the decay constant is given by [9] 

afAB = (am A + am B ) \/WsAq/2 , (2.7) 

where Vs is the spatial volume. Note that bare values of masses and amplitude can be used without 
renormalization thanks to the partially conserved axial current (PCAC) relation of HISQ. Hence, 
the uncertainty from renormalization is absent. 



3. Least-squares Fits 

We use least-squares minimization to find the best fits. The objective function to be minimized 
is the augmented % 2 function, which includes the Gaussian priors as well as usual % 2 , 

xi g (0)= X 2 (0) + ^ { -^-^, (3.1) 

a °a 
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Table 1: Parameters of HISQ ensembles are tabulated. With a few intentional exceptions, the light quark 
masses (m/) are approximately 0.2, 0.1, and 0.037 times the physical strange quark mass. The exceptions 
appear in the ensembles with unphysical values of the strange quark mass and non-degenerate up and down 
quark masses. The HISQ ensembles also include three ensembles to study finite-volume effects, where all 
parameters are the same except the lattice volume. 





light-light 




heavy-light 


No. of states 


J+K = 1+0 




J+K = 2+1 


Priors 


No 


Mi 

K 


-M = 700 ±140 MeV 
-M = 400±200MeV 




<2 (fhl) ^min ^max 




<2 (fhl) ^min ^max 




0.15 16 23 




0.15 8 23 


Fit range 


0.12 20 31 




0.12 10 31 



0.09 30 47 0.09 15 47 

0.06 40 71 0.06 20 71 



Table 2: Our fit methods are summarized. J + K stands for the number of states included in the fits, where 
J (K) is that of ordinary (alternate) states. For heavy-light, the priors are applied to the mass gaps. The 
minimum distances are about 2.4 fm for light-light and 1.2 fm for heavy-light. 
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Figure 1: A example of the fractional errors for various valence quark mass combinations are shown. The 
circle and square symbols stand for random and Coulomb wall sources respectively. Here, "light" stands for 
up and down quarks. 



where 6 represents a set of parameters to be estimated, and a runs over parameters constrained 
with central value pL a and width o a . 

Table 2 summarizes our fit methods. We divide the correlators into two groups depending 
on the valence quark masses: light-light and heavy-light, where by "heavy" we mean charm and 
"light", up, down, and strange. For light-light, we use 1+0 state with large minimum distances 
(~ 2.4 fm). For heavy-light, however, we use 2+1 state with constrained masses and the reduced 
minimum distances (~ 1.2 fm). We choose minimum distances small enough to provide adequate 
statistical signal for included states but not so small as to be polluted from ignored excited states. 

The main reason for the differences in the two groups of fits is the fractional errors (noise- 
to-signal) of the correlators [10]. Figure 1 shows an example of the fractional errors for various 
valence quark combinations. One can easily see that the fractional error grows much faster for 
heavy-light. This means that good signals are not expected at large distances for heavy-light. 
Thus, it is better to reduce the minimum distance to get better statistical precision, while including 
additional heavier states in the fitting function. 

The central values of the priors are taken to be M £ ,( 2 550) ~~ ~ 700 MeV for the ordinary 
excited state and Mjy —Mp « 400 MeV for the alternate state [6]. The widths are chosen narrow 
enough to effectively stabilize the fits, but not so narrow as to significantly influence our determi- 
nation of the ground states, the states of interest. We find that these priors allow us to get good 
confidence levels of fits (or p-values) over all heavy-light fits. 
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Figure 2: The ratios of diagonal elements of the covariance matrix for the random wall source correlator 
on the 0.12 fm physical quark mass ensemble are plotted for various values of block size b. On the left, 
the correlators are averaged over four measurements on each lattice. On the right, however, only a single 
measurement is used to measure the correlators. For details, please see the text. 

4. Autocorrelation of Correlators 

We use the blocking method to compute the correlators and their statistical errors from the 
Monte Carlo samples. The blocking method divides measurements in a time series into contiguous 
blocks and uses their means as independent measurements in successive analysis. As a proxy for 
autocorrelation, we examine ratios of the covariance matrix of the correlators as a function of the 
block size, 

$=CJW/Cg ] . (4.1) 

where b denotes the size of blocks, and cfj represents the covariance matrix computed with the 
block size b. 

Figure 2 shows examples of the diagonal elements of the ratio matrix. On the left plot, we can 
see that the autocorrelation has a peculiar (Euclidean) time dependence. This happens because we 
measure the correlators from four source times and average them to get one estimator on each gauge 
configuration. Then, when we proceed to the next gauge configuration, we move the locations of 
sources by a constant amount Nt/& plus a small displacement 8 = ± 1 , which prevents moving back 
and forth to the same places. Note that the positions of peaks coincide with the "moving distance", 
^7/8 + 8, or its multiple. We can compare these with the results of the single measurement (with 
the the moving distance ~ Nt/4), which are plotted on the right in Fig. 2. Hence, we find these 
autocorrelations are induced by the way we place the multiple sources thoughout the lattices. 

To suppress the effects of the autocorrelations, we need to increase the size of the blocks. 
However, as the block size increases we have fewer blocks, and our error estimates get less accu- 
rate. In this trade-off, we find that an optimal block size is four, which gives us enough statistics 
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to estimate the covariance matrices of large dimension while reducing the underestimation of error 
due to the autocorrelations for the ensembles where the number of configurations is about a thou- 
sand. We use the block size two or one depending on the number of configurations when it is less 
than a thousand. 

5. Excited State Effects 

To address the systematic uncertainty due to ignored excited states, we vary the fitting range 
and priors and then test changes of the values of interest. For measurements of the decay constants 
on the 0.09 fm physical sea quark mass ensemble, for example, no statistically significant deviation 
is found beyond statistical precision with any of the following variations: with reduced prior widths 
(~ 50 MeV) for the mass gaps, with shifts of the central values (~ ±100 MeV), with additional 
constraints on the amplitudes, with changes of the minimum distances (~ ±2), and so on. This 
leads us to our conservative estimations of the systematic uncertainty from the excited states in 
Ref. [1]. 
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